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^ ■ Abstract. We report a controllable method for producing mixed two-photon states 

via Spontaneous Parametric Down- Conversion with a two-type-I crystal geometry. 
By using variable polarization rotators (VPRs), one obtains mixed states of various 
purities and degrees of entanglement depending on the parameters of the VPR. The 
generated states are characterized by quantum state tomography. The experimental 
^^ ■ results are found to be in good agreement with the theory. The method can be 

easily implemented for various experiments which require the generation of states with 
controllable degrees of entanglement or mixedness. 
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1. Introduction 

Quantum information is a promising field that utilizes the nonclassical aspects of 
physical systems to perform sophisticated tasks of computations and communications. 
Quantum entanglement plays an important role in the implementation of these tasks. 
Arguably, the most famous entangled states are the Bell states, which are pure 
and maximally entangled. They are used in quantum cryptography [Ij, quantum 
teleportation [2j, and the demonstration of various concepts of quantum mechanics 
[31 Sj. However, apart from the highly entangled pure states, there exists a vast 
uncharted region of state space, where states can be simultaneously mixed and entangled 
[5l [6]. Mixed states are useful in investigations of quantum computing 0, studies of 
the quantum-classical interface [i8j, and decoherence channels [9]. These applications 
motivate an interest in the generation and characterization of mixed states. 

Up to now, several methods have been suggested for generating photonic mixed 
states. In Ref. [lO], Werner states are produced with controlled addition of white 
noise to the Bell states. In Ref. [llj, more sophisticated schemes for producing broad 
classes of states are presented. Such schemes employ an incoherent temporal mixing of 
state amplitudes, several decoherers, or a hybrid technique. However, as commented in 
Ref. [H], such schemes require many additional optical components, and are challenging 
in practice. 

In this paper, we report a controllable method for producing mixed states, which 
requires only few additional optical components added to a Bell-state generation set-up. 
We use variable polarization rotators ( VPRs) , placed in the pump and signal beams of a 
conventional Spontaneous Parametric Down-Conversion (SPDC) set-up with two type-I 
crystals. Varying the parameters of the VPRs enables us to obtain mixed states of 
various purities and degrees of entanglement, which are then characterized by quantum 
state tomography. We remark that VPRs have been used earlier by Gogo et al in studies 
of quantum erasure p[2j without, however, presenting a detailed study and systematic 
characterization of the generated states. 

2. Mixed-state generation with VPR 

The basic principle underlying this method of producing mixed states is to generate 
incoherent mixtures of orthogonal Bell states with controllable weights. Consider a 
SPDC set-up with two type-I crystals with orthogonal axes, pumped by a laser polarized 
at -45° llSj. The produced state is the BeU state |$_) = (|HsHi) - |VsVi))/V2, 
where H and V denote horizontal and vertical polarizations of down-converted photons, 
respectively, and the subscripts s-signal, and i-idler (hereafter omitted for simplicity) 
label the spatial modes. If the pump is polarized at +45°, the orthogonal Bell state 
|$+) = (|HH) + |VV))/v^ is generated. By blending these two Bell states with different 
weights, one obtains the incoherent mixture 

p = a|$_)($_| + (1 - a)|$+)($+| (1) 
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= i(|HH)(HH| + |VV)(VV|) + (« - ^) (|HH)(VV| + |VV)(HH|), 

where a {0 < a < 1) is the weight of the state |$_)($_| in the mixture. 

To generate such mixed states, one switches the polarization of the pump back and 
forth between —45° and +45° and averages over a sufficient interval of time to sample 
both polarization states. Such a switching can be realized by inserting a VPR with 
electrically driven retardance in the pump beam. By applying a square waveform across 
the VPR, one achieves fast flippings between the two orthogonal linear polarizations. 
The parameter a is equal to the duty cycle (DC = r/T, where r is the duration of the 
high voltage, and T is the period of the waveform) of the applied square waveform. 

The purity P = tr{p^} of such states is given by 

/ 1\2 1 
P = i-^) +2- (2) 

When a is zero or one, one of the Bell states is produced, and the purity is one. When 
a is 0.5, the two Bell states are mixed with the same proportion, and the purity has 
the minimum value of 0.5. Thus, changing a enables one to obtain states with different 
purities. 

Different amounts of entanglement can also be obtained by varying the parameter 
a. The tangle T is a measure of quantum-coherence properties of a quantum state 
[m [15]. It has a value of zero for separable states, and one for maximally entangled 
states. For the states given by Eq. (1), the tangle is related to a by 

r = (1 - 2af. (3) 

Generated states can be characterized by polarization correlation analysis. For the 
Bell states, the polarization correlation gives a visibility of 100% when measured in the 
±45° basis. For the states given by Eq. (1), the visibility varies with a according to 

V = \l-2a\ = Vf. (4) 

In particular, when a = 0.5, the generated state is p = |( |HH)(HH| + |VV)(VV|), 
which results in zero visibility of the polarization correlation in ±45° basis. Note that 
for the state given by Eq. (1), V = a/T, however, this relation is not true in general. 

With a single VPR in the pump beam and a set to 0.5, one produces mixed states 
with minimal purity P = 0.5, see Eq. (2). However, for a completely mixed state, 
P = l/d, where d is the dimension of the density matrix (rf = 4 for a two-photon state). 
The completely mixed state can be written as: 

PM = ^(|HH)(HH| + |VV)(VV| + |HV)(HV| + |VH)(VH|) = ^. (5) 

Thus, one has to include in the mixture the other two Bell states |^_) = (|HV) — 
|VH))/y2, and |^+) = (|HV) + |VH))/y2 with equal weights. These two BeU states can 
be generated in the same set-up by inserting a VPR with a half- wave retardance in either 
the signal or idler beam. In this case the VPR transforms the states |$_) — > l^-)? 
and |$+) — > 1^+)- By switching between zero and half- wave retardance of the 
VPR, incoherent mixtures of |$_)($_| and |^_)(^_|, or |$^)($^| and |^_)(^_|, are 
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produced. Hence, by using two VPRs with duty cycles set to 0.5, one placed in the 
pump beam, and the other in the signal or idler beam, one can obtain the state given 
by Eq. (5). More generally, one can obtain mixtures of non-maximally entangled states 
I^arb) = /3|HH) +7I VV) by rotating the polarization of the pump beam to an arbitrary 
angle. 

3. Experimental Set-up 

3.1. State Preparation 

The experimental set-up is shown in Fig. [H Two beta-barium borate (BBO) crystals of 
0.8 mm thickness with orthogonal optic axes are cut at 29.5° for noncollinear frequency- 
degenerate phase matching. The BBOs are pumped by a frequency stabilized continuous 
wave diode laser at a wavelength of 404.6 nm (Ondax, LM series). The first crystal 
produces pairs of vertically polarized photons from the pump's horizontal polarization 
component, while the second crystal produces pairs of horizontally polarized photons 
from the pump's vertical polarization component. With noncollinear SPDC, the down- 
converted photons travel in different directions with a cone opening angle of 4°. Down- 
converted photons of narrow spatial bandwidth are collected into single-mode optical 
fibers (SMF) using aspherical lenses (L) with a focal length of 11mm, placed at a 
distance of 800 mm from the BBOs. The down-converted photons are filtered with 
interference filters (IF) centered at 810 nm, with 10 nm full width at half maximum. 
Polarization rotations in the optical fibers are compensated for by manual polarization 
controllers (PC). 

With a half- wave plate (HWP) in the pump beam at 22.5°, the vertical polarization 
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Figure 1: Experimental set-up. A cw diode laser pumps two type-I BBO crystals with 
orthogonal axes, where the SPDC occurs in the noncollinear frequency-degenerate regime. 
Mixed states are generated by inserting variable polarization rotators (VPRs) in the pump 
and signal beams. QP are quartz plates used to control the phase of the produced states. M 
are mirrors. The SPDC photons are coupled into single-mode fibers (SMF) with lenses (L). PC 
are polarization controllers, IF are interference filters. Quarter- and half- wave plates (QWP, 
HWP) and polarizing beam splitters (PBS) are used for quantum state characterization. Dl-4 
are single photon detectors, whose outputs are processed by a coincidence circuit (&). 
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of the pump is rotated to —45°, and the resulting state is (|HH) — e^^\YY))/^/2^ where 
(j) is the relative phase. This phase can be set to zero or integer multiples of 27r by 
manipulating the ellipticity of the pump with two quartz plates (QP) which are placed 
after the HWP, thus producing |$_) = (|HH) - |VV))/y2. 

Two different types of polarization rotators have been used in the experiment. 
The first one is a photoelastic modulator (HINDS Instrument, I/FS50). It switches 
between zero and half- wave retardance at a frequency of 50 kHz, with a fixed DC of 
0.5. The second one is a liquid crystal retarder (LCR) (Meadowlark, LRC-200). It is 
configured so that zero retardance is obtained at the "low" voltage of l^ = 1.39 V. The 
pump polarization entering the crystal remains at —45°, and |$_) is generated. For the 
"high" voltage of V^ = 1.64 V, half-wave retardance is obtained. The pump polarization 
entering the crystal changes to +45°, and |$+) is generated. A square waveform with 
the frequency of 1 Hz is applied to the LCR. The averaging is done over 3-minutes time 
intervals. By varying the duty cycle of the square waveform, different mixed two-photon 
states are obtained, see Eq. (1). 

To produce the completely mixed state, given by Eq. (5), the photoelastic 
modulator is placed in the pump beam, and the LCR with DC = 0.5 is placed in 
the signal beam. The operation voltages of the LCR for down-converted photons are 
adjusted to Vi, = 1.84 V and Va = 6.95 V. The frequency of the LCR square waveform 
is kept at 1 Hz, and the averaging is done over 3-minutes time intervals. 

3.2. State characterization 

To fully characterize the states, quantum state tomography is performed. In both signal 
and idler arms, we install quarter-, and half- wave plates (QWP, HWP) and a polarizing 
beam splitter (PBS) [I6]. The transmission and reflection ports of both PBSs are 
directed to single-photon detectors (Silicon Avalanche Photodiodes, quantum efficiency 
^ 50%, Qutools Twin QuTD). Coincidence events are registered using the Time- To- 
Digital Converter (quTAU) with a time window of 5ns. The coincidences between any 
two of the detectors are recorded. By manipulating the wave plates in front of the 
PBSs, we perform measurements in nine different bases, which give an overcomplete 
measurement of 36 outcomes. Using the technique of maximum likelihood estimation 
[Uj, the state is inferred from the collected data. 

4. Results 

For the calibration of the set-up, the polarization correlation analysis in the ±45° basis 
is performed with the Bell state |$-). The constant "low" voltage is applied across the 
LCR in the pump beam. The QWPs in the signal and idler arms are fixed at 0°. The 
HWP in the idler arm is fixed at —22.5° and the HWP in the signal arm is rotated. 
This causes the coincidence rate between the counts of two detectors in the transmitted 

iV_|_45o — iV_45o 

ports of the PBSs to vary sinusoidally. The visibility is defined as 1/ = 



A^+45° + A^-45° 
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Figure 2: Dependence of the visibility 
of the polarization correlation measure- 
ments in the ±45° basis on the DC of the 
LCR (solid circle), and of the photoelastic 
modulator (red open diamond). The solid 
line is the theoretical prediction. The er- 
ror bars are smaller than the symbols. 
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Table 1: Fidelities of the reconstructed states with 
the target states for various DCs of the LCR. F > 
97% are consistently obtained for all the states. The 
last two entries with PEM and 1/4 refer to the state 
generated using only the photoelastic modulator in 
the pump beam, and the completely mixed state 
generated using two VPRs, respectively. 



where A^+450 and A^_45o are the coincidence rates when the HWP in the signal arm is 
oriented at +22.5° and —22.5°, respectively. The obtained value of visibility 97.3±1.4% 
indicates high-quality of the produced Bell state. The analysis is then extended to 
the mixed states prepared with the LCR operating at different DCs, and with the 
photoelastic modulator. The obtained experimental visibilities are in good agreement 
with the theoretical expectation of Eq. (3), see Fig. [2l One observes a degradation of 
the visibility as a increases from to 0.5. After which, the visibility is gradually restored 
as a increases from 0.5 to 1. 

Next, quantum state tomography on the generated states is performed. The matrix 
elements of reconstructed density matrices are shown graphically in Fig. 3, for a=0.05, 
0.25, and 0.5. The fidelities, F, of the reconstructed density matrices with the states 
given by Eq. (1) are calculated. Using the definition F = tT{\y/py/a\}^ where p is 
the maximum likelihood estimator, and a is the target state, fidelities F > 97% are 
consistently observed, see Table. [H 

As an example, for a = 0.25, the reconstructed density matrix expressed in the HV 

basis is 

0.545 0.049 + 0.0121 

0.049 - 0.012i 0.008 
-0.013 - 0.038i -0.005 - O.OOSi 
-0.232 -O.llOi 0.015-0.0041 
The reconstructed state has a tangle of T = 0.2476 ± 0.0049, a purity of P = 
0.6295 ± 0.0022, and a fidelity with the target state of F = 0.9814 ± 0.0006. 

As one can see from Fig. 3, an increase of a from to 0.5 causes vanishing of 
the off-diagonal elements, which indicates the decreasing amount of entanglement of 
the mixed states. The calculated values of the tangle for the reconstructed states, are 
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Figure 3: Absolute values of real (a, b, c) and imaginary (d, e, f) parts of the density matrices 
representing the reconstructed states for 0.05DC (a, d); 0.25DC (b, e) and 0.50DC (c, f). The 
vanishing of the off-diagonal elements is clearly seen. 
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Figure 4: Dependence of tangle (a) and purity (b) of reconstructed states on the duty cycle 
of the LCR (solid circles), and of the photoelastic modulator (red diamond). The solid curves 
are the theoretical predictions. The error bars are smaller than the symbols used for both 
figures. 



in good agreement with Eq. (3), see Fig. [l](a). The state with almost zero tangle is 
generated with the photoelastic modulator. It performs better than the LCR at a=0.5, 
as the wavelength of the pump laser is at the edge of the working range of the LCR. 
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Figure 5: Absolute values of real (a) and imaginary (b) parts of the reconstructed density 
matrix for the completely mixed state. 



Lastly, one verifies that the method indeed produces states with different amounts 
of mixedness. The dependence of the purity of the states on a is shown in Fig. [11(b). 
The obtained results are close to those predicted by Eq. (2). 

To demonstrate generation of the completely mixed state, the set-up is modified by 
inserting VPRs in the pump and in the signal beams, see Sec. 3.1. The reconstructed 
density matrix of the state is shown graphically in Fig. [5l The reconstructed state has 
a tangle of zero, a purity of P = 0.2615 ± 0.0007, and a fidelity with the target state of 
F = 0.9942 ± 0.0003. Such a state is separable and highly mixed. 

5. Conclusion 

In summary, an accessible way of generating controllable mixtures of the Bell states with 
high fidelities is reported. There is no need for extra optical components, except the 
VPRs. Furthermore, the method does not require any reconfiguration of the standard 
set-up for Bell states generation, as the VPRs can be easily turned on and off, to generate 
mixed and pure states, respectively. In addition, it offers the important advantage of 
interferometric stability. The presented method will facilitate state generation and can 
be used in various quantum-information processing applications. We are planning to 
use it for an experimental realization of qubit-pair tomography with witness bases [18j . 
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